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Abstract 



In this paper we study general nonlinear stochastic differential equations, where the usual 
Brownian motion is replaced by a Levy process. We also suppose that the coefficient mul- 
tiplying the increments of this process is merely Lipschitz continuous and not necessarily 
Ph linear in the time-marginals of the solution as is the case in the classical McKean-Vlasov 

model. We first study existence, uniqueness and particle approximations for these stochastic 
differential equations. When the driving process is a pure jump Levy process with a smooth 
but unbounded Levy measure, we develop a stochastic calculus of variations to prove that the 
time-marginals of the solutions are absolutely continuous with respect to the Lebesgue mea- 
sure. In the case of a symmetric stable driving process, we deduce the existence of a function 
solution to a nonlinear integro-differential equation involving the fractional Laplacian. 

Key words: Particle systems; Propagation of chaos; Nonlinear stochastic differential 
equations driven by Levy processes; Partial differential equation with fractional Laplacian; 
Porous medium equation; McKean-Vlasov model. 
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<N 

C ' This paper studies the following nonlinear stochastic differential equation: 

O " 

O : {x t = Xo + J t a(X s -,P s )dZ s , t€[0,T\, 

1 Vs G [0, T], P s denotes the probability distribution of X s . 

• i— i 

We assume that Xq is a random variable with values in IR fe , distributed according to m, (Z t )t<T 
a Levy process with values in R d , independent of Xq, and a : R k x V(R k ) -» R kxd , where V(R k ) 
denotes the set of probability measures on R k . Notice that the classical McKean-Vlasov model, 
studied for instance in [22], is obtained as a special case of ([1]) by choosing a linear in the second 
variable and Z% = (t,Bt), with E>t being a (d — l)-dimensional standard Brownian motion. 

The first section of the paper is devoted to the existence problem and particle approxima- 
tions for dH). Initially, we address the case of square integrable both, the initial condition Xq, 
and the Levy process (Z t )t<T- Under these assumptions the existence and uniqueness problem 
for dU) can be handled exactly as in the Brownian case Zt = (t,Bt). The nonlinear stochas- 
tic differential equation ([T]) admits a unique solution as soon as a is Lipschitz continuous on 
R k xV2(R k ) endowed with the product of the canonical metric on R k and the Vaserstein metric d 
on the set V2(R k ) of probability measures with finite second order moments. This assumption is 
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much weaker than the assumptions imposed on a in the classical McKean-Vlasov model, where 
it is also supposed to be linear in its second variable, that is, a(x,v) = L fe s(x,y)u(dy), for a 
Lipschitz continuous function ? : M fc x M. k — > R fcxrf . Then, replacing the nonlinearity by the 
related interaction, we define systems of n interacting particles. In the limit n — > +00, we prove, 
by a trajectorial propagation of chaos result, that the dynamics of each particle approximates the 
one given by (pQ). Unlike in the very specific McKean-Vlasov model, where the universal C/^/n 
rate of convergence corresponds to the central limit theorem, under our general assumptions on 
a, the rate of convergence turns out to depend on the spatial dimension k. 

In the next step, the square integrability assumption is relaxed. However, to compensate 
for its loss, we assume a reinforced Lischitz continuity of a : the Vaserstein metric d on V(M. k ) 
is replaced by its smaller and bounded modification d\ defined below. Then, choosing square 
integrable approximants of the initial variable and the Levy process, we prove existence for ([I]). 
Uniqueness remains an open question. 

In the second section, we deal with the issue of absolute continuity of Pt when Z is a 
pure jump Levy process with infinite intensity. For the sake of simplicity of the exposition, 
we restrict ourselves to the one-dimensional case k = d = 1. When a does not vanish and 
admits two bounded derivatives with respect to its first variable, and the Levy measure of Z 
satisfies some technical conditions, we prove that, for each t > 0, Pt has a density with respect 
to the Lebesgue measure on M. The proof depends on a stochastic calculus of variations for 
the SDEs driven by Z which we develop by generalizing the approach of Bichteler-Jacod [4], 
(see also Bismut [5j), who dealt with the case of homogeneous processes with a jump measure 
equal to the Lebesgue measure. In our case, the nonlinearity induces an inhomogeneity in time 
and the jump measure is much more general, which introduces additional difficulties making 
the extension nontrivial. Graham-Meleard [TT] developed similar techniques for a very specific 
stochastic differential equation related to the Kac equation. In that case, the jumps of the process 
were bounded. In our case, unbounded jumps are allowed and we deal with the resulting possible 
lack of integrability of the process X by an appropriate conditioning. 

In the third section, we keep the assumptions made on a in the second section, and 
assume that the driving Levy process Z is symmetric and a-stable. Then, we apply the absolute 
continuity results obtained in Section 2 to prove that the solutions to ([I]) are such that for 
t > 0, Pt admits a density pt with respect to the Lebesgue measure on the real line. In addition, 
calculating explicitly the adjoint of the generator of X, we conclude that the function pt(x) is a 
weak solution to the nonlinear Fokker-Planck equation 

d tPt (x) = D«(\a(., Pt )\ a Pt (.))(x) 
lim t ^ + p t (x)dx = m(dx), 

where, by a slight abuse of the notation, a(.,pt) stands for a(.,pt(y)dy), the limit is understood 
in the sense of the narrow convergence, and = — (— A) Q//2 denotes the spatial, spherically 
symmetric fractional derivative of order a defined here as a singular integral operator, 

D«/(x) =K I (f(x + y) - f(x) - l M<1] f{x)y) 

Jr \y\ 

where K is a positive constant. For 

1 _^ 

a(x, v) = [g e * v{x)) s with e > 0, g £ {x) = e 2 = and s > 0, 

V27T£ 

one obtains the nonlocal approximation dtPt = Dx((9e*Pt) as Pt) of the fractional porous medium 
equation dtPt = D^(p" s+1 ), the physical interest of which is discussed at the end of the paper. 



2 



Other nonlinear evolution equations involving generators of Levy processes, such as fractional 
conservation laws have been studied via probabilistic tools in, e.g., [13] . and[H]. 

Notations : Throughout the paper, C will denote a constant which may change from 
line to line. In spaces with finite dimension, the Euclidian norm is denoted by | |. Let V(R k ) 
denote the set of probability measures on R k , and V 2 (R k ) ~ the subset of measures with finite 
second order moments. For /i, v G V2(R k ), the Vaserstein metric is defined by the formula, 

d(n, v) = inf { ( f \x- y\ 2 Q(dx, dy)"] : Q G V(R k x R k ) with marginals /x and v 1 . 



It induces the topology of weak convergence together with convergence of moments up to order 
2. The modified Vaserstein metric on V(R k ) defined by the formula, 



v) = inf {[ / \x- y\ 2 A 1 Q{dx, dy)) : Q 6 T(R k x R k ) with marginals fi and v 

[ \JR k xR k J 

simply induces the topology of weak convergence. 



1 Existence of a nonlinear process 

We first address the case when both, the initial condition Xq and Z are square integrable, before 
relaxing these integrability conditions later on. 



1.1 The square integrable case 



In this subsection we assume that the initial condition Xq, and the Levy process (Zt)t<T, are 
both square integrable : E(|Xo| 2 + |.Zt| 2 ) < +oo. Under this assumption, the following inequality 
generalizes the Brownian case (see, |20| . Theorem 66, p. 339) : 



Lemma 1 Let p > 2 be such that E(\Zt\ p ) < +oo. There is a constant C p such that, for any 
R kxd -valued process {Ht)t<T predictable for the filtration Tt = o~(Xq, (Z s ) s <t), V£ £ [0,7"], 



E ( sup 



H u dZ u 



<C P [' E(\H s \ p )ds. 
Jo 



Because of this inequality for p = 2, the results obtained for the classical McKean-Vlasov 
model driven by a standard Brownian motion still hold. First, we state and prove 



Proposition 2 Assume that Xq and {Zt)t<T are square integrable, and that the mapping a is 
Lipschitz continuous when endowed with the product of the canonical topology 

on R k and the Vaserstein metric d on V2{R k )- Then equation ([1]) admits a unique solution 
such that E (sup 4<T \Xt\ ) < +oo. Moreover, if for some p > 2, E(|Xo| p + \Zt\ p ) < +oo, then 
E (sup t<T \X t \ p ) < +oo. 



Proof : We generalize here the pathwise fixed point approach well known in the classical 
McKean-Vlasov case (see, Sznitman [22]). Let B denote the space of cadlag functions from [0, T] 



3 



to R k , and 7> 2 (B) the space of probability measures Q on D such that J B sup t<T \ Y t \ 2 Q(dY) < 
+00. Endowed with the Vaser stein metric Dt(P, Q) where, for t <T, 



D t (P, Q) = inf \ (( sup \Y S - W s \ 2 R(dY, dW)) : R G V(B x D) with marginals P and Q 1 , 

[ \JDxID) s<t J J 

T^O^) is a complete space. 

For Q G "P(B) with time-marginals (Qt)tg[o,T] > m view of Lebesgue's Theorem, the distance 

d(Qt,Q s )< [ \Y t -Y s \ 2 Q{dY) 
Jo 

converges to 0, as s decreases to t (respectively, d(Q t -,Q s ) < f n \Y t - — Y s \ 2 Q(dY) converges to 
0, as s increases to t; here Q t - = Q o l^Z 1 is the weak limit of Q s as s —> t~). Therefore, the 
mapping t G [0, T] — > Q t is cadlag when T^O^) is endowed with the metric d. As a consequence, 
for fixed x G the mapping t G [0, T] — > a(x,Qt) is cadlag. Hence, by a multidimensional 
version of Theorem 6, p. 249, in [20J, the standard stochastic differential equation 



rt 

X? = X + a{X® , Q s )dZ s , t G [0, T] 
Jo 



admits a unique solution. 

Let denote the mapping on ^(D) which associates the law of X® with Q. Let us check 
that $ takes its values in Vi(p). For K > 0, we set = inf{s < T : \X?\ >K}. By Lemma Q] 
and the Lipschitz property of a, one has 



E (sup|X s Q ATA .| 2 ) < C (e(\X \ 2 ) + J E (l {s < TK} |a(X s Q ,Q s ) - a(0,5 )| 2 + k(0,5 )| 2 ) ^ 

<C(E(|X | 2 ) + / e(sup\X^ Atk \ 2 ) ds + t [ sup\Y t \ 2 Q(dY) + t\a(0,5 )\ 2 ) 
y Jo \r<s J Jnt<T J 

By Gronwall's Lemma, one deduces that 

E (sup|X s Q ArA .| 2 ] <C[E(\X \ 2 ) + \a(0,5 )\ 2 + [ sup \Y t \ 2 Q(dY) ) , 
ys<T y y Jot<T J 

where the constant C does not depend on K. Letting K tend to +00, one concludes by Fatou's 
Lemma that 

f sup|y s | 2 c^(Q)(y)=E(sup|Af | 2 ) < C [E(\X \ 2 ) + \a(0,5 )\ 2 + [ sup \Y t \ 2 Q(dY) ) . 

JHs<T \s<T J y Jot<T J 

(2) 

Observe that a process (-Xt)tg[o,T]> such that E (sup t<T |A" 4 | 2 ) < +00, solves (HJ if and only if its 
law is a fixed-point of <3?. So, to complete the proof of the Proposition, it suffices to check that 
admits a unique fixed point. 

By a formal computation, which can be made rigorous by a localization procedure similar 
to the one utilized above, for P,Q G T^O^) one has 



E 



sup\X? - X?\ 2 ) <C ! E(\a(X^,P s )-a(Xf_,Q s )\ 2 )ds 

s<t J J0 

< C I E (sup |X r p - X®\ 2 \ + d 2 (P s , Q s )ds. 

Jo \r<s J 
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By Gronwall's Lemma, one deduces that, Vi < T, 



E ( sup \Xf - XfA <C f d 2 (P s , Q s )ds. 

\s<t J JO 



Since D t 2 ($(P), $(Q)) < E(sup s < t |Xf - X?| 2 ), and d(P s ,Q s ) < D S (P,Q), the last inequality 
implies, Vi < T, 



Df($(P),$(Q)) < C f D 2 s (P,Q)ds. 
Jo 



'0 



Iterating this inequality, and denoting by the iV-fold composition of we obtain that, 
ViV G N*, 



D*($ N (P),* N (Q)) < c N j Q D 2 s (P,Q)ds < ^-dUp,Q). 

Hence, for iV large enough, & N is a contraction which entails that admits a unique fixed point. 

If, for some p > 2, E(|Xo| p + |^t| p ) < +00, a reasoning similar to the one used in the 
derivation of (j2|, easily leads to the conclusion that the constructed solution {X t )t<T of ([I]) is 
such that 

E fsup \X s \p\ < C ( E(|Xo| p ) + |<7(0, 5 )r + E ('sup l^] 2 ^ | < +00. 



Our next step is to study pathwise particle approximations for the nonlinear process. Let 
((Xq, Z*))j G fsj* denote a sequence of independent pairs with (Xq, Z l ) distributed like (Xq, Z). For 
each i > 1, let (-2Q)te[o,T] denote the solution given by Proposition [2] of the nonlinear stochastic 
differential equation starting from Xq and driven by Z l : 



X l + f t a(X^,P s )dZl, te[0,T] 
Vs G [0, T], P s denotes the probability distribution of X* 



(3) 



Replacing the nonlinearity by interaction, we introduce the following system of n interacting 
particles 



Xi + jla{X^,n n s _)dZl te[0,T], l<i<n, 
where ji n = ^ Y^j=i &xi> n denotes the empirical measure 



(4) 



Since for £ = (x%, 



), and C = (2/1, •■ -,y n ) in 



ink 



one has 



n n 



< 



1 n 

n ^ 

j'=i 



1/2 



n 



ie — ci- 



(5) 



Existence of a unique solution to (j3|), with finite second order moments, follows from Theorem 
7, p. 253, in [20]. Our next result establishes the trajectorial propagation of chaos result for the 
interacting particle system 
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Theorem 3 Under the assumptions of Proposition 

lim supE [ sup \Xl' n -Xf\ 2 ) =0 

Moreover, under additional assumptions, the following two explicit estimates hold : 

• //E(|X | fc+5 + \Z T \ k+5 ) < +oo, then 

supE ( sup \X l t ' n - X\\ 2 J < Cn~T^; (6) 

i<n \t<T J 

• If <t(x, v) = J Kfe y)v(dy), where ? : M fc x IR fc — > M fcxci zs a Lipschitz continuous function, 
then 

S upE(sup|Xr-X t f ] <-, (7) 

where the constant C does not depend on n. 
The proof of the first assertion relies on the following 

Lemma 4 Let v be a probability measure on M fc such that f Rk \x\ 2 v(dx) < +oo , and v n = 
— X]?=i ^tj denote the empirical measure associated with a sequence (£i)i>i of independent ran- 



1 

n £-13 -- 

dom variables with law v. Then, Vn > 1, 



E(c2 2 (z/>)) < 4 f \x\ 2 v{dx), and lim E id 2 (u n , v)) = 0. 

J Rk n-*+oo 



Proof of Lemma [4] : By the strong law of large numbers, as n tends to oo, almost surely 
v n converges weakly to v and, Vi, j € {1, . . . , k}, J Rk XiV n {dx) (resp. j Rk XiXjis n (dx)) converges 
to f Rk Xiiy(dx) ( resp. XiXji'(dx)) . Since the Vaserstein distance d induces the topology of 
simultaneous weak convergence and convergence of moments up to order 2, one deduces that 
almost surely, d(u n ,u) converges to 0, as n tends to oo. Hence, to conclude the proof of the 
first assertion, it is enough to check that the random variables (d 2 (u n , v)) n >i are uniformly 
integrable. To see that note the inequality 



d 2 



The right-hand side is nonnegative and converges almost surely to 4 f Rk \x\ 2 v(dx), as n — > oo. 
Since its expectation is constant, and equal to the expectation of the limit, one deduces that the 
convergence is also in L . As a consequence, for n > 1, the random variables in the right-hand 
side, and therefore in the left-hand side, are uniformly integrable. I 
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Proof of Theorem [3] : Let P n = \ YJj 

=1 5 X j denote the empirical measure of the independent 
nonlinear processes ([3]). By a formal computation, which can be made rigorous by a localization 
argument similar to the one made in the proof of Proposition [21 one has, Vi < T, 

¥,Lu V \XT-XlA <C f\{\a{Xr^ n s )-a{Xl,P^)\ 2 )ds 

\s<t J JO 



+ C f\{\a(XlP^)-a(XlP s )\ 2 )ds. 
Jo 



In view of the Lipschitz property of <r, the estimate ([5]), and the exchangeability of the couples 
(X 1 , X hn )i<i< n , the first term of the right is smaller than C J Q *E ^sup r<s \Xp n — X % r \ 2 ^j ds. By 
GronwalFs Lemma, and the Lipschitz assumption on <r, one deduces that 

E ^sup|X*' n -X\\ 2 ^j < C J\(\a {XlP?) -a{Xi,P s )\ 2 ^j ds < C j\{d 2 (P™ , P s ))ds. 

The first assertion then follows from Lemma [H the upper-bounds of the second order moments 
given in Proposition and by Lebesgue's Theorem. 

The second assertion is deduced from the upper-bounds for moments of order k + 5 com- 
bined with the following restatement of Theorem 10.2.6 in [21] : 



e {d 2 (pj\ p.)) < c (i + )J f Rk \y\ k+5 P* 



(dy) J n k + 4 . 



where the constant C only depends on k. The precise dependence of the upper-bound on 
Jijfc \y\ k+5 Ps(dy) comes from a carefull reading of the proof given in [21j . 

Finally, if, as in the usual McKean-Vlasov framework (see [22]), a(x, v) = f Rk y)v(dy), 
where? = (s a b)a<k,b<d :l*xM* -» R kxd is Lipschitz continuous, then E (|<r (X*, P") - <r(X*, P s )| 2 ) 
is equal to 

k d 



a=l 6=1 j,l=l 



S ab (Xi,Xi)- / <; ab (Xi,y)P s (dy) 



S ab (XiX l s )- / s ab (Xi,y)P s (dy) 



Since, by independence of the random variables Xj , . . . , X" with common law P s , the expecta- 
tion in the above summation vanishes as soon as j ^ I, the third assertion of Theorem 3 easily 
follows. 1 



Remark 5 

• Observe the lower estimate 

d(v n ,v)> [ / min \£j — x\ 2 v(dx) ) > inf ( / min \y~ — x\ 2 v(dx) ) 
\Vi<j<« J (yi,-,yn)em n \JRki<j<n J 

Moreover, according to the Bucklew and Wise Theorem [Bj, if v has a density <p with 

respect to the Lebesgue measure on M k which belongs to L^+ k (M k ), then, as n tends to 
infinity, 

1/2 



n l / k m f / / mm \y. _ x\ 2 ip(x)dx ) 
(yi,-,yn)m k ) n \JR k l <i< n J 
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converges to Cfelly?]! h , where the constant only depends on k. Hence, one cannot 

2 + k 

expect E(<i 2 \v n , u)) to vanish quicker than Cn~ 2 l k . 

Therefore, if v — > cr(x, v) is merely Lipschitz continuous for the Vaserstein metric, one 
cannot expect E ^sup t <^ \X l t ' n — X\\ 2 ^j to vanish quicker than Cn~ 2 l k . The rate of con- 
vergence obtained in ([6]) is not far from being optimal at least for a large spatial dimension 
k. Nevertheless, in the McKean-Vlasov framework, where the structure of a is very spe- 
cific, one can overcome this dependence of the convergence rate on the dimension k, and 
recover the usual central limit theorem rate. 

• The square integrability assumption on the initial variable X$ can be relaxed if a is Lip- 
schitz continuous with M. k x P(R fc ) endowed with the product of the canonical topology 
on M. k and the modified Vaserstein metric d\ on V(M. k ). Indeed, one may then adapt the 
fixed-point approach in the proof of Proposition [2] by defining V as the space of probability 
measures on B, and replacing sup s<t \ Y S — W s \ 2 by sup s<t \Y S — W s \ 2 A 1 in the definition of 
Dt(P, Q). This way, one obtains that the nonlinear stochastic differential equation ([1]) still 
admits a unique solution even if the initial condition Xq is not square integrable. More- 
over, for any probability measure v on R fc , if v n is defined as above, K(d 2 (u n ,u)) remains 
bounded by one, and converges to as n tends to infinity. Therefore, the first assertion in 
Theorem [3] still holds. 

The next subsection is devoted to the more complicated case when the square integrability 
assumption on the Levy process {Z t )t<T is also relaxed. 

1.2 The general case 

In this section, we impose no integrability conditions, either on the initial condition Xq, or on the 
Levy process {Z t )t<T- Let us denote by m G V(M. ) the distribution of the former. According 
to the Levy-Khintchine formula, the infinitesimal generator of the latter can be written, for 
/ G C b 2 (R d ), in the form 




where a = {a>ij)\<ij<d is a non-negative symmetric matrix, b a given vector in M, d , and a 
measure on M. d satisfying the integrability condition / Kd (l A \y\ 2 )f3{dy) < +oo. 

To deal with non square integrable sources of randomness, we impose a stronger continuity 
condition on a, namely, we assume that a is Lipschitz continuous when R k x V(R k ) is endowed 
with the product of the canonical topology on M. k and the modified Vaserstein metric d\ on 
V(M> k ). Notice that, under this assumption, for each x G M. k , the mapping v G V(M. k ) — > a(x, v) 
is bounded. 

In order to prove existence of a weak solution to (pQ), we introduce a cutoff parameter 
N G N*, and define a square integrable initial random variable Xq = Xol^ Xo \<N}i an d a square 
integrable Levy process (Z^) t <T by removing the jumps of (Z t )t<T larger than N : 

z t = z t ~ 1{|AZ S |>7V}AZ S . 

s<t 
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Let (^t V )te[o,r] denote the solution given by Proposition [5] of the nonlinear stochastic differential 
equation starting from Xq and driven by (^/ V )ie[o,T] : 

(x t N = X» + jj a(X^ , P s N )dZ», t G [0, T] ^ 
1 Vs G [0, T], Pg denotes the probability distribution of X^ 

We are going to prove that when the cutoff parameter N tends to do, then (^f V )tg[o,T] converges 
in law to a weak solution of ([I]). More precisely, let us denote by P N the distribution of 
(X t Ar ) tg r 0] r|, and by (^)te[o,Tl the canonical process on D. 

Proposition 6 The set of probability measures (P N )^eN* is tight when D is endowed with 
the Skorohod topology. In addition, any weak limit P, with time marginals (Pt)te[o,T]> °f its 
converging subsequences solves the following martingale problem : 



Pq = m and \/(p : M. k — > R, C 2 with compact support, 
Mf = (p(Yt) — <p(Yo) — Jq C[P s ]ip(Y s )ds J is a P -martingale 

v / L t~~ U.J 



(9) 



where for each v G P(R fc ) ; and any x G M. k , 
1 k 

C[v]<p(x) =- i (Ja(J *( x ^))ij d l l ,x J i P{ x ) + (cr{x,u)b).V(p(x) 

+ / [f{x + a{x,u)y) -if(x) - l {M < 1} a(x,u)y.V(p(x)] (3(dy)ds. (10) 



Proof : Let us first remark that for iV G N*, and for a fixed xGl*, the mapping t G [0, T] — > 
<r(x, P/') is cadlag and bounded by a constant not depending on N. As a consequence, according 
to Theorem 6 p. 249 |20| . for a fixed M G N*, the stochastic differential equation 



X N,M = X NAM + /* a{X ^ M ,P s N )dZ^ AM , t G [0, T], 

Jo 

admits a unique solution. Let us denote by P N > M the law of (X t ' M )te[o,T\- By trajectorial 
uniqueness, ViV G N*, Vt G [0, T], 

V N _ V N,M 
A t — A t > 

as long as |Xo| Vsup fe ] 0) T] |AZ(| < M. The probability of the latter event tends to one as M tends 
to infinity. Using both the necessary and the sufficient conditions of Prokhorov's Theorem, one 
deduces that the tightness of the sequence (P^jjvgff* is implied by the tightness of the sequence 
(P N > M ) Nm *, for any fixed M G N*. 

Let us now prove this last result by fixing M G N*. Using the boundedness of d\, one 
easily checks that 

sup E ( sup \Xt' M \ 2 ) < +oo. (11) 

N£N* \t<T J 

This implies tightness of the laws of the random variables \swp t<T \X^' M \) . In order to 

use Aldous' criterion, we set e, S > 0, and introduce two stopping times S, and S, such that 
< S < S < (S+5)AT. Let us also remark that, for K G N*, and b K = b+ f Rd yl {1< \ y \< K yf3(dy), 
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the process ( Zf = Zf - b K t 



observe that 



te[o,T] 



is a centered Levy process and therefore a martingale. Now, 



xf M -Xf M \ 2 >e\ <I 



+ : 



l S a(X?' M ,P s N )b NAM ds 
Js 



> 



a(X N > M ,P?)dZ? AM 



e 



(12) 



Using the boundedness of the sequence (& )iveN*j the Lipschitz property of a with respect 
to its first variable, and (jllh combined with the inequalities of Markov and Cauchy-Schwarz, 
one obtains that the first term of the right-hand-side is smaller than C5 2 /e, where the constant 
C does not depend on N . For the second term of the right-hand-side, one remarks that Doob's 
optional sampling Theorem, followed by the Lipschitz property of a, and (jlip . imply that 



E 



a{X N,M : pN )dZ; 



NAM 



E 



(T 0,(7 



N,M r)N\ 



+ 



N,M 

S J 



P. 



y\ 2 i{\ y \<NAM}P(dy) 



ds 



< C6, 



where C does not depend on N. By Markov's Inequality, the second term of the right-hand-side 
of (|12p is smaller than CS/e and, in view of Aldous' criterion, we conclude that the sequence 
(P N ' M ) NeN * is tight. 

Now, let us denote by P the limit of a converging subsequence of (P n )n£N* that we still 
index by iV for simplicity's sake. Also, let ip denote a compactly supported C 2 function on R k . 
For p £ N*, < s\ < S2 < ■ ■ ■ < s p < s < t < T, and a continuous and bounded function 
ip : (R k ) p — ► R , let F denote the mapping on V{p) defined by 



F(Q) = J L(Y t ) - <p(Y s ) - £ £[Q u ](p(Y u )dv) 4>(Y S1 ,. . .,Y Sp )Q(dY). 

For F N defined like F, but with ln y i <N j(3(dy) replacing j3{dy) in the definition ([TO]) of C[v], one 
has F N (P N ) = 0. Therefore 



\F(P 



F{P N )-F 1V (P 1 ")\ <2(t-s 



N t t>N\ 



J 



l{\y\>N}(3{dy) 



N- 



0. 



The mapping (x,u) G R x "P(R ) — > C[u]ip(x) is bounded, continuous in x for a fixed v, and 
continuous in is, uniformly for x G R fc . Therefore, as soon as si, . . . , s p , s,t do not belong to the 
at most countable set {u €]0, T] : P(AY U / 0) > 0}, then F is continuous and bounded at point 
P which implies F(P) = limjv^ +00 F(P N ) = 0. In view of the right continuity of u — > Y u and 
Lebesgue's theorem, this equality still holds without any restriction on s±,... ,s p ,s,t. Hence, 
(M^) tg [ 0> T] is a -P-martingale. Since the sequence (Xq)^^* converges in distribution to Xq, 
Pq = m, which concludes the proof. I 



The above existence result for the martingale problem implies an analogous existence 
statement for the corresponding nonlinear Fokker-Planck equation. 
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Proposition 7 Let P denote a solution of ([9]). Then the time marginals (Pt)t£[o,T] solve the 
initial value problem 

d t P t = C*[P t ]P u P = m, (13) 

in the weak sense, where, for v G V(^. k ), C*[u] denotes the formal adjoint of C[u] defined by the 
following condition: V</>, ip C 2 with compact support on M. k , 

/ £* [v]ip(x)(p(x)dx = / ip{x)C[v\<p{x)dx. 

Moreover, the standard stochastic differential equation 

Xf = X + f a(X^,P s )dZ s (14) 
Jo 

admits a unique solution (^t P )tG[o,T] an d, for each t £ [0, T], Xf is distributed according to Pt. 

Proof : The first assertion follows readily from the constancy of the expectation of the P- 
martingale (Af^)tG[o,T] • Existence and uniqueness for (|14|) follows from [20], Theorem 6 p. 249. 
Now, if Q t denotes the law of Xf for t > 0, then (Qt)t>o solves 

d t Q t = C*[P t ]Qt, Qo = m : 

in the weak sense. Since (Pt)t>o also is a weak solution of this linear equation, by Theorem 5.2 
[3], one concludes that Vi < T, P t = Qt- 1 



Remark 8 We have not been able to prove uniqueness for the nonlinear martingale problem 
Q. However, our assumptions, and Theorem 6, p. 249, in [20J, ensure existence and uniqueness 
for the particle system Like in the proof of Proposition [6l one can check that the laws of the 
processes X ,n , n > 1, are tight. According to [22j . this implies uniform tightness of the laws 
7r n of the empirical measures /i ra . For a fixed x G M. k , the function v — > o(x,v} is continuous 
and bounded when ^(R^) is endowed with the weak convergence topology. Then one can prove 
that the limit points of the sequence (vr n ) n give full weight to the solutions of the nonlinear 
martingale problem Q. 

2 Absolute continuity of the marginals 

In this section we restrict ourselves to the one-dimensional case k = d = 1, and assume that 
Z is a pure jump Levy process with a Levy measure j3 which admits a density, say j3\, in the 
neighborhood of the origin, that is 

P(dy) = i^KlPi^dy + l| y |>i/3(dy). 
We set 0i(y) = 0, for \y\ > 1. Then 

Z t = ( yN 1 (ds,dy)+ [ y N 2 (ds,dy), (15) 

J(0,t]xM J(0,t]xR 

where Ni, and iVg, are two independent Poisson point measures on R + x M with intensity 
measures equal, respectively, to \\ y \ < i(5i{y)dy, and l\ y \ > iP(dy), and N± is the compensated 
martingale measure of N\. 

We work here 
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• either under the assumptions of Subsection 1.1, i.e., E(|Xo| 2 + \Zt\ 2 ) < +00 and Lipschitz 
continuity of a on R X V2 (R) endowed with the product of the Euclidean metric and the 
Vaserstein metric, 

• or with the general assumptions of Subsection 1.2, i.e., no integrability conditions on Z and 
Xq, and Lipschitz continuity of a on R x 'P(R) endowed with the product of the Euclidean 
metric and the modified Vaserstein metric d\. 

Propositions [2] and [6] ensure the existence of a probability measure solution P of ([9|) . According 
to Proposition [7J there is a unique pathwise solution X (which is then unique in law) , to the 
stochastic differential equation 

X t = X + [ a(X s -,P s )yN 1 (dy,ds)+ [ a(X s - , P s ) yN 2 (dy, ds), (16) 
J(0,t]xR J(0,t]xR 

and, for t £ [0,T], Xt is distributed according to the time marginal Pt. 

Roughly speaking, our goal is to prove that, for each t > 0, the probability measure Pt has 
a density with respect to the Lebesgue measure as long as the measure f3, restricted to [—1,1], 
has an infinite total mass due to an explosion of the density function (5\{y) at 0. Indeed, we 
have in this case an accumulation of small jumps immediately after time 0, which will imply the 
absolute continuity of the law of Xt under suitable regularity assumptions on /3\ . 

For this purpose we develop a stochastic calculus of variations for diffusions with jumps 
driven by the Levy process defined in (fT5j) . We thus generalize the approach developed by 
Bichteler and Jacod [I] (also, see Bismut [5]), who considered homogeneous processes with a jump 
measure equal to the Lebesgue measure. Here, the nonlinearity introduces an inhomogeneity in 
time, and the jump measure is much more general, which complicates the situation considerably 
and introduces additional difficulties. Graham and Meleard [11] developed similar techniques 
for a very specific stochastic differential equation related to the Kac equation. In that case, the 
jumps of the process were bounded. In our case, unbounded jumps are also allowed. 

Our approach requires that we make the following standing assumptions on the coefficient 
a(x,v), and the Levy density f3\\ 

Hypotheses (H): 

1. The coefficient a(x,f) is twice differentiable in x. 

2. There exist constants K\, and K 2 , such that, for each x, and v, 

\<j' x {x, v)\ < Ki, and \a"(x,u)\ < K 2 . (17) 

3. For each x, and is, 

a(x,u)^0. (18) 

Hypotheses (Hi): 

1. The function (3\ is twice continuously differentiable away from {0}. 
2. 

J Pi(v)dy = +00. (19) 

3. There exists a non-negative and non-constant function k of class C 1 on [—1, 1] such 
that k(— 1) = k(l) = 0, and such that 
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J 1 k\y)p l {y)dy 



< +00, 



(20) 



for all y £ [-1,1], 



sup — | 



P 1 (y + X(l + ay)k(y)) 



My) 



(l + \(ak(y) + (l + ay)k'(y))-l 



1 



(21) 



sup 



(1 + ay)k{y) + ak{y) + (1 + ay)k'{y) 



Pi(y)dy < +00, 



(22) 



sup 

ae[-K 



f 1 ( ^(y + A(l + ay)%)) , 
P / sup -r-^ (1 + A(a%) + (1 + ay)k (y)) 

Ae[-i,i] V Puy) 



+ 



(3[(y + X(l + ay)k(y)) 



k 2 {y)f5 1 {y)dy < +00, 



k 2 (y) 



(23) 



• for all y £ [-1, 1], 

|kte)l< i(rTifo ; * iirbo- (24) 

The assumption (Hi 3) on /?i is obviously technical, but the assumption (Hi 2) is essential, 
and cannot be avoided if one hopes to prove the absolute continuity result. The main example 
satisfying assumptions (Hi) is the symmetric stable process with index a £ (0,2), for which 
0{dy) = K dy/\y\ 1+a , as developed in Section 3. 



Theorem 9 Consider the real-valued process X satisfying the nonlinear stochastic differential 
equation (I16p . Assume that a satisfies Hypotheses (H), and that j3\ satisfies Hypotheses (Hi). 
Then the law of the real-valued random variable X? has a density with respect to the Lebesgue 
measure. 



The remainder of this section is devoted to the proof of Theorem [9] which will proceed 
through a series of lemmas and propositions. Our aim is to show that Pt has a density with 
respect to the Lebesgue's measure. Because of the compensated martingale term N\ it would 
be natural to work with square integrable processes. But the fmiteness of E (sup t<r |Xj| 2 ) is 
not guaranteed because of the big jumps of ./V2 . So, to develop the relevant stochastic calculus 
of variations in L 2 , we use a trick defining Ft as the conditional law of (Xq, Ni, N%) given 
(Xo,iV^), where N2 denotes the restriction of the measure N2 to [0, T] x R. Thus, in what 
follows, the random variables considered being functions of (Xq, N±, N2) we may define their 
Ey-expectations as the integral of the corresponding functions under Ft- From now on, for 
notational simplicity, every statement concerning Ey, or Ft, holds almost everywhere under the 
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law of (X , A 2 T ), even if this fact is not mentioned explicitly. This conditioning allows us to use 
the same techniques as if the process X were square integrable. More precisely, given Aj, there 
are finitely many jump times and jump amplitudes of N2 on (0, T] and we will denote them by 
(Xij Y{), ■ • • , (Tfc, Y k ). 



Lemma 10 

E T (sup|A t | 2 ) < +00. (25) 

t<T 



Proof : As usual, we localize via r n = inf{t > 0, |AQ| > n}. Then, for t <T, 

E T (sup|A tAr J 2 ) < c[\X \ 2 + J ^ \y\ 2 (E T (\X s \ 2 )+a 2 (0,P s ^j(3x(y)dyds 

k . 
+ £ \Y\ 2 (Et(|X Ti _| 2 ) + a 2 (0, P T J) l Ti < tATn 
i=l ' 

|X | 2 + sup a 2 (0,P u )( 1 + V|li| 2 ) + / " E T (\X s \ 2 )ds 
ue[o,T] \ ~[ J h 

, rt/\T„ r \ 

+ sup|y i | 2 / E T (|X S _| 2 ) / N 2 (dy,ds) . 
i=l JO ^|y|>l / 

At this point we apply Gronwall's Lemma in its generalized form (with respect to the measure 
ds + f^yx ^(dy, ds), see, for example, Ethier and Kurz |10j . p. 498). The result then follows. 



Let us now explain our strategy to prove Theorem [9] before giving the technical details. We 
are going to prove that there exists an a.s. positive random variable DXt such that Mt(DXt) < 
+00, and G C~(R), 

\&T{<t>'{X T )DX T )\ < CUWoo, (26) 

where C£°(R) denotes the space of infinitely differentiable functions with compact support on 
the real line. Indeed, in the special case DXt = 1 this inequality implies that the conditional 
law of Xt given (Xq, Aj) admits a density (see, for example, Nualart [17], p. 79). If DXt 7^ 1, 
the law of Xt under Qt = eJ(DXt) a£ imits a density. But since DXt > 0, a.s., then is 
equivalent to Pt, and the conditional law of Xt given (Xq, N£) still admits a density. Of course, 
this implies that the law Pt of Xt admits a density. 

We will prove inequality f|26[) employing the stochastic calculus of variations. Consider 
perturbed paths of the process on the time interval [0, T] and introduce a parameter A 6 [—1, 1], 
sufficiently close to 0. The perturbed Poisson measure N± will satisfy N® = Aq, and be such 
that, for a well chosen P-p-martingale {G^)t<T, the law of its restriction to [0, T] under G t .¥t 
will be equal to the law of Aq under Pt. The process X x will be defined like X, only replacing 
Nx by N± in the stochastic differential equation. Then, for sufficiently smooth functions <f>, we 
will have 

E T (0(X T )) = Et(Gt0(*£)). ( 27 ) 
Differentiating in A at A = 0, in a sense that we yet have to define, we will obtain 

Et{<P'{X t )DXt) = -E T (DG T (t>(X T )), (28) 



14 



where DXt = -^X t \\ = q, an d DGt = -§\G t \\=q. Then one easily deduces (f2"6"j) with C = 
E r (|/X; r |). 

Let us describe the perturbation we are interested in. Let g be an increasing function of 
class C£°(]R), equal to x on [— i; |], equal to 1, for x > 1, and to —1, for a; < —1. Note that 
Hslloo < 1; an d that g(x)x > 0, for x 6 R*. 

Now, we define the predictable function v : f2 x [0,T] x [—1, 1] i— > K via the formula 

«(*,!/) = l {s>S} (l + y^(X s _,P s )) g(a(X s ^,P s )) k(y) (29) 

where 5 is a stopping time that we are going to choose later on in order to ensure that DXt > 0, 
a.s. (see the discussion before Proposition [T5]) . It is easy to verify that the function y \— > v(t,y) 
is of class C 1 on [—1, 1], and in what follows we will denote its derivative by v'(t,y), Also, for 
every u, t, and y, 

\v(t,y)\ < (l + KtMv) and \v'(t, y)\ < K x k{y) + (1 + Ki)\X(y)\. (30) 

For A G [—1,1], let us introduce the perturbation function 

7 A (t,y) =y + Xv(t,y). (31) 

We can easily check that, for every u>, and t, the map y \— > j x (t,y) is an increasing bijection 
from [-1, 1] into itself, since by and ([30]), \v'\ < ~, and k(-l) = k(l) = 0. 

Let us denote by N x the image measure of the Poisson point measure N± via the mapping 
7 A defined, for any Borel subset A of [0,T] x [—1, 1], by the integral 



N?(A) = J lAfa-ffayVN^dt). 
We also introduce the function 

#i(v) 

which appears below in the definition of the process G x (in Proposition fT2|) . As a preliminary 
step, we obtain the following estimates concerning V x . 

Lemma 11 There exists a constant C such that, for almost all uj, and for all s £ [0, T], 

sup \\V\s,y)-l\ < \ , (33) 

A,3/G[-l,l] A 1 



i r 1 

sup — / \v x (s, y )-i\ 2 My)dy <c, 

= [-1,11 A 



Ae[- 

1 

sup - 4 , 

Ae[-i,i] A J-i 



V x (s,y)-l-\-^V x (s,y)/ x= o 2 



(34) 



/3i(y)% < C. (35) 



Proof : Inequality (j33H follows from (|21j) . Also, one has 



(36) 
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and 



(s.W = v (s,y){l + \v {s,y)) -^-j- + 2 — v(s,y)v (s,y). 



Since, for A E [—1, 1], we have estimates, 



V\s,y)-l 



< 2 



d 

dX 



V\s,y)/x =0 



+ 



A 2 



V x (s,y)-l-X—V x ( S ,y)/ x=0 



and 



V x (s,y)-1-X— V x (s,y)/ X=0 



2 A 4 

< — sup 
4 



one deduces d35j) (resp. ([MI)) from ([23]), and PH)(resp. (1221). ([23]) . and (121)1 



(37) 



We are now ready to introduce the promised earlier definition of the process G x . 
Proposition 12 (%) For every A £ [—1,1], i/ie stochastic differential equation 
G X = l+ [ GUV X (s,y)-l)N 1 (dy,ds), 

J(0,t]xR 

has a unique solution G x which is a strictly positive martingale under Ft and such that 



(38) 



sup E T I sup|G A | 2 ) < +oo. 
Xel-i,i] \t<T J 

(ii) The law of under F T = G T .F T is the same as the law of N\ 



(39) 



under '. 



Proof : (i) Thanks to ((34]), the stochastic integral M 4 A = J (Q t]xK (V rA (s, y) - l)N x (dy, ds) is well 
defined and is a Fx square integrable martingale. The unique solution to ([38]) is the exponential 
martingale G t A = S{M x ) t = e M * rio<s<t( 1 + ^)e" 4M ' given by the Doleans-Dade formula. 

Using (|33|) . we remark that the jumps of M x are more than —1/2 so that G x is positive 
for each t. Moreover, using (|34p . Doob's inequality and Gronwall's Lemma, we deduce from (|38[) 
that (131 holds. 



(ii) Let us denote fx(dy, dt) = li y i <l /3i(y)dydt, and compute the image measure 7 (V ./x). 
For a Borel subset A of [0,T] x [—1, 1], we have 



l A (t, y + Xv(t, y))V x {t, y)(3 l (y)dydt 
l A (t,y')^l(l + Xv'(t,ym(y)- * l{ 



Hv) v ' a,n LKa 'l + \v'{t,y) 

l A {fy')(5 1 {y')dy'dt = fi(A), 



dt 



(40) 



where y' = y + Xv(t, y). Hence 



1 X {V X .lf)=LL. 



(41) 
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Since N± is independent from (Xq, N%), the compensator of N± under P^ is /x. By the Girsanov's 
theorem for random measures (cf. Jacod and Shiryaev [12], p. 157), its compensator under 
p£ = Gy.Pr is V x .fi and thus, the compensator of N x = 7 A (iVi) is equal to 7 A (F A ./i) = fi. 
We have thus proven that the compensator of N x under Pj, is fj,, and the second assertion in 
the proposition follows. I 



Next, we study the differentiability of G x with respect to the parameter A, at A = 0. 
Proposition 13 (i) The process 

DG t = [ ± V \s,y) /x=0 N 1 (dy,ds)= [ (v'(a,y) + ^\v(s,y)) N^ds) 

(42) 

is well defined, and such that 

E T (sup|L>G t | 2 ) < +oo. (43) 

t<T 

(ii) The process G x is L 2 -differentiable at A = ; with the derivative DG which is under- 
stood in the following sense: 



E T sup \G X - 1 - \DG t \ 2 = o(A 2 ), a.s., (44) 
\t<T J 

as A tends to 0. 



Proof : (i) Thanks to ((36]) and ([22]), for almost all to, and all s G [0,T], 



A 2 



Tx V\s,y) /x 



Pi(v)dy < C. (45) 



'-l 

Therefore the process DGt is well defined and satisfies ([33 
(ii) Moreover, one has 



E T sup|G A - 1 - XDG t \ 2 
\t<T 



< C 



< C 



Er ((g x (V x (s,y) - l) - ^V x (s, y) /A=0 ) j f3 1 (y)dyds 
^E T ^G x (v x (s,y)-l-X-^V x (s,y) /x=0 ^ j Px{y)dyds 
J (ot]xR E T ^(^ yA ( s ^)/A=o)) (G X S ~ l) 2 ^ Hv)dyds. 



+\ 2 C 



Now, according to (|35p and (|39p . we obtain that 



^E T UG x (V x (s,y)-l-X-^V x (s,y) /x=Q ^j \ f) x (y)dyds < C\H. 
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Furthermore, by (Ho]) . 

l(0t]xR ET ((^^^'^^O ~ 1)2 ) " i^ 7 ^ " 1)2)& ' 

and using ([34"]) and ([39]) . we may show that, for each t < T, 

E T ((G X - l) 2 ) = / E T ({G x (V x (s,y) - l)f)f3 x {y)dyds < CX 2 . 

J(0,t]xR 



This concludes the proof. 



In the next step we define the perturbed stochastic differential equation. Let us recall that 
the probability measures P± are fixed and are considered as time-dependent parameters. Thus 
the process X is a function Fp(Xq, Ni, N 2 ) of the triplet (Xq, Ni, N 2 ). 

Define X x := F P (X , N? t N 2 ). Hence, using Proposition [12] (ii), the law of X x under 
is equal to the law of X under Py. A simple computation shows that X x is a solution of the 
stochastic differential equation 



X? = X + 



y a(X x _,P s ) (Nf(dy, ds) - Pi{dy)ds) + f y a(X x _, P s )N 2 (dy, ds) 

(0,(]xl v 7 J(0,t]xR 

X + / (y + \v(s,y)) a(X x _,P s ) (x^dy^s) - V x (s,y)(3 1 (dy)ds) 

I In +1 ssi® \ / 



+ / y a(X x _,P s )N 2 (d yj ds), ( since 7 A (^V) = A*) 

J(0,t]xR 

= X + [ ya(X x _,P s )N 1 (dy,ds) + X [ a(X x _, P s )v(s, y)Ni(dy, ds) 

J(0,t]xR J(0,(]xR 

+ f ya(X x _,P s )N 2 (dy,ds)- [ (y + Xv(s,y))a(X x _, P s )(V x (s,y) - \)^{y)dyd8. 

J(0,t]xR J(0,t]xM 

(46) 

Using ([20]) for the second term, the fact that j\{y 2 + k 2 (y))(5\(y)dy < +00, ([34"]) . and the 
Cauchy-Schwarz inequality for the last term, we easily prove that equation (|46p has a unique 
pathwise solution. 

Let us now show that X x is differentiable in A, at A = 0, in the L 2 -sense. More precisely 
we have the following 



Proposition 14 



(0 Pr ysup \X t x -X t \ 4 j < CX A . (47) 
(ii) P T ( sup \Xj" -X t - XDX t \ 2 ) = o(X 2 ), (48) 

\t<T 



as A tends to 0, where DX is a solution of the affine stochastic differential equation 

DX t = [ ya' x (X s -,P s )DX s „ N 1 (dy,ds) + I a(X s _, P a )v(a,y)Ni(dy, ds) 
JtOMxR J(0,t]xR 
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+ / ya' x (X s -,P s )DX s - N 2 (dy,ds)- f y a{X s -,P s )(p x {y)v' {s,y) + P[(y)v(s,y))dyds. 

J(0,t]xM J(0,t]x[-l,l] 



(49) 



Proof : In order to prove assertion (i), we need the following moment estimate : 



sup E T sup|X t A | 4 < +00. 
Ae[-i,i] \t<T J 

It relies on the following classical estimation (see, [3J, Lemme (A. 14)) 

„2/ \ o f.,u„, \ 1 / 4 



(50) 



Er I I / H s p(y)N 1 (dy,ds) ) \<C 

'(0,i]xR / 



Ay)/h(v)dy) + J ^(y)f3 l {y)dy 



X Et I sup |if u | ) ds, 

/o \u<s 



(51) 



for any predictable process H, and any measurable function p : [—1,1] > R such that the right- 
hand side is finite. Conditioning by iVj', the times and the amplitudes of jumps of N2 on (0, t] 
are given by Yi), • • • , (T fc , and 

A: 



E 



'T 



ytr(X2_,P 8 )N 2 (dy,ds) 



<C^F 4 4 (E T (|^_| 4 ) + f T 4 (0,P r j) 



<Csup|F| 4 [ sup a 4 (0,P n )+ / / E T (|X s A _| 4 )iV2(rfy,^) . 
«=i \«e[o,T] J J|y|>i / 

Applying (|5ip with /?(y) = y, and p(y) = fc(y), (|34p and Gronwall's Lemma with respect to the 
measure ds + f i >:l N 2 (dy, ds), we easily check (|50j) and deduce 

sup E T f sup|a(X t A , J P t )| 4 J < +00. 
Ae[-i,i] \*- T / 

Now, we write X x — Xt using (|16p and (|46p . Assertion (%) is obtained following an analogous 
argument. 

To prove (ii) we need to isolate the term Z x = X x — Xt — XDXt , and as in [4] , Theorem 
(A. 10), we write 

X? -X t - XDX t 

y Z^o^X^PsXmidy, ds) + N 2 (dy, ds)) 

(0,i]xM 

+ f y U(X^,P S ) - a(X s .,P s ) - a' x (X s „,P s )(X x _ - X.-j) (Ni(dy, ds) + N 2 (dy, ds)) 
+ / Xv(s,y) (a(X x _,P s )-a(X s ^,P s )) Nx(dy,ds) 

J(0,t]xR v 7 

y a(X s ,P s ) (v x (s,y) - 1 - X^-V x (s, y) /x=0 ) PMdyds 

(0,i]xR V aA J 

y (a(X x ,P s )-a(X s ,P s )) (V x (s,y) - l)^(y)dyds 



(0,t]xR 



'(0,t]xR 



X v(s,y)a(X x ,P s )(V x (s,y) - l)^{y)dyds. 
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Under Hypotheses (H) and (Hi) all the integral terms, except the first one, are of order 
A 2 . Indeed, for the second term, we use Taylor's expansion of a, and ([37]) : for the third term, 
we use (j20j) . and ((171) ; for the fourth we use the Cauchy-Schwarz inequality, and ((351) ; for the 
fifth term, we use the Cauchy-Schwarz inequality, (|34p . and (|47|) : for the sixth term, we use 
Cauchy-Schwarz inequality, (f34l) . and (f20|) . Then, as previously, using Gronwall's Lemma for 
the conditional expectation, we obtain the result. I 



The term DXt requires our special attention. Observe that, after integration by parts (in 
the variable y), the last term in (I14|) cancels the compensated part of 

/ CT(X s _,P>(s,y)iVi(dy,ds), 

J(0,t]xR 



and one obtains 



DX t = I DX s _dK s + L t (52) 
J o 



K t = [ y a' x (X s _,P s )N 1 (dy,ds) + I y <(X S _, P s ) N 2 (dy , ds) , (53) 

J(0,t]xK ./(0,t]xR 



where 



= / a(X s _,P s )v(s,y) Nxidy^s). (54) 

J(0,t]xR 

As in [12J, Theorem 4.61, p. 59, or in [3], we can solve (|52p explicitely. The jumps Aif s are of 
the form ya' x (X s , P s ). Thus 1 + AK S may be equal to and then the Doleans-Dade exponential 

£{K) t = e Kt U 0<s < t (l + AK s )e~ AK ° 

vanishes from the first time when AK S = — 1. We follow [4] to show that DXt ^ 0, but the 
strict positivity (which has not been proved in the latter) necessitates a careful analysis. 

Let us define the sequence of stopping times 5i = inf{t > 0, AK t < —1}, Sk = inf{i > 
Sk-i,AKt < —1}, So = 0. Since a' x is bounded, there is a finite number of big jumps on the 
time interval [0, T], so that there exists an n such that S n <T < S n+ i = do, and P(S n = T) = 0. 

Solving equation (f52j) gives 

DX t = £(K - K Sk ) t [DX Sk + [ (1 + AKs^SiK - K Sk );ldL s ) 

V J(s k ,t] J 

if Sk < t < •S'fe+i, and t <T, (55) 
where ^ 5fc = Ks k w In particular, 

DXt = £(K - K s ") T (DX Sn + [ (1 + AK S )-^(K - ^ 5 ");_ 1 dL s ') . 

V V(S„,T] / 

Because of the definition of SVi, the exponential martingale £{K — K Sn ) s is non-negative on 
[S n ,T]. If the perturbation v did not vanish before time S n , it would not be clear how to 
control the sign of DX$„ ■ That is why we choose S = S n in ([29]) : 



v(s,y) = l {s>Sn} (l+ya' x (X s ^,P s )) k(y) g(a(X s ^,P s )) 
so that DX$ n = 0. For this choice, we obtain 
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Proposition 15 We have DXt > 0, almost surely. 

Proof : One has DX T = £(K - K s ") T Y n , where Y n = f.^ T] (l + AK S )~ 1 £(K - K s ™)-}dL s . 

Since Ni and N% are independent, the sets of jumps are almost surely distinct and then 
1 + AK S can be replaced by 1 every time the jump of K comes from a jump of iVjj. So, 



Let us consider the set A n = {(u, s, y), h n (u, s, y) > 0} and define the stopping time r = ini{t > 
S nJ(s n ,t]xl-i,l] h n (s,y)Ni(ds,dy) > 0} = inf{t > S n , J {Snit]x[ _ hl] lA n (s,y)Ni(ds,dy) > 0}. 

Using the definitions of v and 5 n +i, one knows that, if S n (u>) < s < T A 5 n+ i(a;), then 

(u,s,y) G A n 44> o-(X s _(uj),P s ) / 0, 
which is always the case in view of Hypothesis (H). 

On the other hand, /( gn)T ] x [_ 1 ^ ^■A n (s,y)N 1 (ds,dy) < 1, 
so K(J, S _ 1 1 | 1a„ (s, y)Ni(ds, dy) J < 1 and, for a; in a set of probability 1, 



These two remarks, and the fact the J_ 1 f3\{y)dy = +oo, imply that r = S n , almost surely. So 
Y n is strictly positive. Therefore DXt is strictly positive as well and the proof is complete. I 



We are now in a position to complete the proof of Theorem [9l 

Proof of Theorem 10 : For 4> £ C£°(M), we differentiate the expression (f2T|) in the L 2 -sense 
with respect to A, at A = 0, and hence, we obtain the "integration-by-parts" formula (f25j) . Then, 
since E,t\DGt\ < +oo, we obtain (|26p . which concludes the proof. I 

3 The case of a symmetric stable driving process Z 

In this section, we assume that the real- valued driving process Z is a symmetric stable process 
with index a £ (0, 2), i.e., Z is given by ([To]) with f3{dy) = rmra dy, where K > is a normaliza- 
tion constant. The generator of this process is the fractional Laplacian (or, fractional symmetric 
derivative) of order aonl: 




where 



K(s,y) 



£(K - K S ")~}(1 + ya f x (X s ^,P s )y 1 v(s,y)o-(X s „,P s 
£(K - K s -);}k(y)a(X s .,P s )g(a(X s ^,P s )) > 0. 




££/(*) = K I (f(x + y) - f(x) - l m < 1} f(x)y) 



dy 




l + Q ' 



This operator may be defined alternatively via the Fourier transform T : 
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When <t:Ix V(M) — > R satisfies Hypotheses (H), it is possible to explicitly calculate the adjoint 
C*\y\ involved in the nonlinear Fokker-Planck equation (|13|) . For smooth functions (p,tp ■ R i— > R, 
one has 

£[i/]<p(a;)i/'(a:)da; = iT / [ip{x + s(x)y) - ip(x) - l { \ y \< 1} s{x)ytp' \xj) ip{x)dx, 
Jr 2 ' \y\ 

where s(x) = a(x,u). Setting z = -s(x)y, and observing that f R {l{\ z \< s (x)} ~ 1{| 2 |<i}) = 
0, one gets 

dy 



(<p(x + s{x)y) - <p(x) - l { \ y i< 1} s(x)yip'(x) \ 



dz 

(<p(x -z)- ip(x) + l { | z |<i } V(^)) \ s (. x )\ a -TJT+^- 

J IK \ Z \ 

Since 

/ (p(x — z)[\s\ a ifj](x)dx = / ip(x)[\s\ a '4'](x + z)dx, 
Jr ' Jw 

and 

/ (p'(x)[\s\ a ^](x)dx = - [ ip(x)[\s\ a ij]'(x)dx, 
Jr Jr 

invoking Fubini's theorem one concludes that 

/ C[v]<p(x)ip(x)dx = K [ <p(x) [ f[| S rV](x + z)-[| S rV](x)-l {k |< 1} z[| S rV] / (x)')^^dx. 
Jr Jr Jr \ / \ z \ 

Therefore 

C*[u}^(x)=D-(\a(.,u)\^(.))(x). 

Moreover, the absolute continuity result given in Theorem [9] permits us to prove existence of a 
function solution to the nonlinear Fokker-Planck equation. 

Theorem 16 Let m G V(M.), and a G (0,2). Assume that the function o~{x,v) satisfies hy- 
potheses (H) and is Lipschitz continuous in its second variable when V(M) is endowed with the 
modified Vaserstein metric d%. Then, there exists a function (t,x) G (0, T] x R i— ► pt(x) G R+ 
such that, for each t G (0, T], j R pt(x)dx = 1 and, in the weak sense, 

(d tPt (x) = D%(\o-(.,p t Tp t (.))(x) (5g) 

I lim t _ >0 + pt{x)dx = m(dx) for the weak convergence, 
where, by a slight abuse of notation, a(.,pt) stands for a (.,pt(y)dy). 



Proof : Existence of a measure solution {Pt)te[o,T] to the nonlinear Fokker-Planck equation 
follows from Propositions[6]and[71 So to conclude the proof, it is enough to exhibit a perturbation 

Then, by Theorem [21 for 



K 



L {M<i}M^ 



function k{y) satisfying hypotheses (Hi) with /3i(y) 
each t G (0, T], we have Pt = pt{x)dx . 

For 7 > 2f, and e G (0, 1/2), let k e denote the even function on [—1, 1] defined by 



k £ (y) 



,1+7 



£ l+7 + (1 + 7 ) £ 7(y _ e ) _ (1 + c )(„ _ £ ) 

{ (1 + 7 - c)e 1+ ^ - c(l + 7 )e 7 (z/ - 2e), 



1+7 



for y G [0,e], 
for y G [e, 2s] 
for y G [2e,l], 
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where c 



(i+^)-e(i+2-y) ' so * n& t ke(l) = 0. The function k e is non-negative and C 1 on [0,1], 
satisfies {20]) and, My G [-1, 1], 



My) < ( 2 + 7)^ 1+7 , WM\ < (1 + 7) max(l, c)e^, and 



fc e (y) 



<(l + 7)er (57) 



In particular, for small enough e, fc e satisfes (|21|) . Since 



fl(v). 



/?i(y) 



1 + ay)k £ {y) + afc £ (y) + (1 + ay)k' £ (y) 



Pi(y)<c 



k 2 e (y) 



+ k 2 (y) + (k' £ ) 2 (y) 



~C'\y\-^ a -^\ 

in the neighbourhood of 0, (|22p is satisfied as well. In the same way, in the neighbourhood of 0, 

I3'{{y + \{l + ay)k e {y)) l 



sup 

ae[-K u Ki],\€[-l,l] 



+ 



13'^y + \{l + ay)k £ {y)) 



-(1 + \(ak £ (y) + (1 + ay)k' £ (y)) 



My) 



k 2 MMy)<c 



|2+2 7 I 



y 4 y 2 



and (|23j) is satisfied. 

Finally, for a G [-Ki,^] and y,\ G [-1,1], by ([57]), for e < ((1 + «i)(l + 7))" 1/7 , 



1 



/9i(y + A(l + oy)A; e (2/)) 



/?i(y) 



(l + A(afc 6 (y) + (l + aj/)*£(y))-l 



1 |1- |l + A(l + ay) 



k s (y) tl+a 

y 



+ X(ak £ (y) + (1 + ay)k'M)\ 



|A| 



|l + A( l + 02/) feM|i+c 



< 



1 



(1-(1 + Ki)(l+ 7 )^) 1+Q 
and (p^|) is also satisfied for small enough e 



(1 + a )(l + (1 + #i)(l + 7)e 7 )°(l + ^i)(l + 7)^ 7 

+ Kx{2 + 7 )e 1+7 + (1 + #i)(l + 7) max(l, c)e 7 



Remark 17 One of the motivations for our work was to generalize the probabilistic approxi- 
mation of the porous medium equation 

d tPt {x) = D 2 x (p q t (x)), q > 1, (58) 

developed, among others, by Jourdain [9] to the fractional case where D 2 is replaced by D". 
The equation (|58p . which describes percolation of gases through porous media, and which is 
usually derived by combining the power type equation of state relating pressure to gas density 
p, conservation of mass law, and so called Darcy's law describing the local gas velocity as the 
gradient of pressure, goes back, at least, to the 1930's (see, e.g., Muskat [IS]). The major steps 
in the development of the mathematical theory of (|58p were the discovery of the family of its 
self-similar solutions by Barenblatt (see, pQ, and [2]) who obtained this equation in the context 
of heat propagation at the initial stages of a nuclear explosion, and an elegant uniqueness result 
for (|58p proved by Brezis and Crandall [7J. A summary of some of the newer developments in 
the area of the standard porous medium equation can be found in a survey by Otto |18j . 
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However, in a number of recent physical papers, an argument was made that some of 
the fractional scaling observed in flows-in-porous-media phenomena cannot be modeled in the 
framework of (I58p . In particular, Meerschaert, Benson and Baeumer [15 J replace the Lapla- 
cian D%. in (|58p by the fractional Laplacian D" while considering the linear case (q = 1) in a 
multidimensional case of anomalous (mostly geophysical) diffusion in porous media, while Park, 
Kleinfelter and Cushman [19] continue in this tradition and derive scaling laws and (linear) 
Fokker-Planck equations for 3-dimensional porous media with fractal mesoscale. 

On the other hand, Tsallis and Bukman [23] suggest an alternative approach to the anoma- 
lous scaling problem (in porous media, surface growth, and certain biological phenomena) and 
consider an equation of the general form 

d t p\{x) = -D x {F{x)f t (x)) + Dliplix)), r,q£R, (59) 

where F(x) is an external force. The authors manage to find exact solutions for this class of 
equations using ingeniously the concept of Renyi (-Tsallis) entropy but, significantly, suggest 
in the conclusion of their paper that it would be desirable to develop physically significant 
models for which further unification can possibly be achieved by considering the generic case of 
a nonlinear Fokker-Planck-like equation with fractional derivatives. This is what we endeavored 
to do taking as our criterion of " physicality" the existence of an approximating interacting 
particle scheme. For the most obvious, simply-minded generalization, dtPt(x) = D"(p^(x)), that 
physical interpretation seems to be missing, or, at least, we were unable to produce it and, as a 
result, our study lead us to settle on an equation like ([56]) . Indeed for 

1 _£_ 

a(x, v) = (g e * v{x)) s with e > 0, g e (x) = e ^ and s > 0, 

\/27T£ 

(j56j) writes dtpt = D°((g £ *pt) as pt) which, for now, we are viewing as a "physically justifiable", 
fractional, and strongly nonlinear "extension" of the classical porous medium equation. Of 
course, this is only the beginning of the effort to understand these types of models. 
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